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Universal Non-adiabatic Holonomic Gates in Quantum Dots and Single-Molecule 
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Geometric manipulation of a quantum system offers a method for fast, universal, and robust 
quantum information processing. Here, we propose a scheme for universal all-geometric quantum 
computation using non-adiabatic quantum holonomies. We propose three different realizations of 
the scheme based on an unconventional use of quantum dot and single-molecule magnet devices, 
which offer promising scalability and robust efficiency. 
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Holonomic quantum computation (HQC), originally 
conceived by Zanardi and Rasetti [l|, has become one 
of the key approaches to perform robust quantum com- 
putation. The idea of HQC is based on the Wilczek-Zee 
holonomy [2[ that generalizes the Berry phase [3| to non- 
Abelian (non-commuting) geometric phases accompany- 
ing adiabatic evolution. It was shown in Ref. [l[ that 
adiabatic quantum holonomies generically allows for uni- 
versal quantum computation. Different physical settings 
for adiabatic HQC have been proposed 14]. Recently, uni- 
versal HQC based on Anandan's work [5( on non-Abelian 
generalization of the non-adiabatic geometric phase [6[ 
has been proposed [7| . This scheme allows for high-speed 
implementations of quantum gates; a feature that can 
be particularly useful in solid-state devices, where qubits 
typically have short coherence time. Non-adiabatic HQC 
requires fewer energy level transitions than in the corre- 
sponding adiabatic settings 7]. Simulations have shown 
that non-adiabatic HQC can be made robust to a wide 
class of errors, such as decay, dephasing, and fluctuations 
in control parameters [8J|. 

In this Letter, we develop a scheme for conditional 
holonomic gates in a four-level configuration. We illus- 
trate how non-adiabatic holonomic gates in our scheme 
can be implemented in solid state devices consisting of 
few-electron quantum dots |9| and single-molecule mag- 
nets (SSMs) [lOj, which are promising candidates for 
quantum computation [111 ] . By utilizing the scalability 
of these systems, the geometric nature of the proposed 
holonomic gates is a feature that may allow for robust 
and fast manipulations of a large number of qubits. 

The general structure of our system is described by a 
four-dimensional effective state space whose dynamics is 
governed by a Hamiltonian of the form 



H(t) = m(t) i T ° \ 



, T = 



T\\ T12 
T21 T22 



(1) 



and satisfies det T ^ 0, where the latter ensures that T 
is invertible. The Hamiltonian is turned on and off by 
means of the time-dependent scaling function f2(t). We 
demonstrate in the following how H(t) can be used for 
universal non-adiabatic HQC and how it can be realized 
in systems of four coupled quantum dots or in triangular 
molecular antiferromagnets. 

The Hamiltonian in Eq. ([T]) naturally splits the state 
space into two subspaces spanned by A4q = {\a} , |6}} and 
■Mi = {\ c ) i\d)}- Non-adiabatic holonomy transforma- 
tions are realized if (i) the Hamiltonian vanishes on each 
evolving subspace and (ii) the projection operators Pq = 
\a) (a\ + \b) (b\ and Pi = \c) (c\ + \d) (d\ move in a cyclic 
fashion. Condition (i) is satisfied since Pi(t)H(t)Pi(t) = 
U(t,0)PiH(t)PiW(t,0) = 0, I = 0,1, where the first 
equality follows from the fact that H (t) commutes with 
the time evolution operator U(t,0) = e~^/ h ^a H ( s ) ds . 
For condition (ii), we first note that since T is assumed 
to be invertible, there is a unique singular value decom- 
position T = UqDUI, where D = diag{a, (3} and Uo,U\ 
are the positive and unitary parts, respectively. By us- 
ing the singular value decomposition, the time evolution 
operator can be expressed as 



U(t,0) 



_ ( U cos (a t D) Ul -iU sin (a t D) u\ 
-iUi sin (a t D) U x cos (a t D) u\ 



(2) 



in the ordered orthonormal basis M. = {\a) , \b) , |c) , \d)}. 
We assume that T is complex-valued, time-independent, 



where a t = J Q ft(s)ds. The evolution of the P;'s is thus 
cyclic provided there is a t — r such that sin (a T D) = 
diag{sin(a r a), sin(a r /3)} = 0. Under this condition the 
holonomy matrices for the two subspaces I = 0, 1 read 

U{d) = UiZ n u}, n = 0,l, (3) 

where cos (a T D) = Z n with Z = diagjl, — 1}. Co 
and Ci are the paths of the two subspaces spanned by 
M.q and M.\, respectively. These paths reside in the 
Grassmannian manifold G(4;2), being the space of all 
two-dimensional subspaces of the four-dimensional state 
space. One can see from Eq. ([3]) that the holonomy is 
non-trivial if and only if n = 1. 
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Next, we encode the basis states in M. into the com- 
putational two-qubit states \a) — > 1 00) , \b) — > |01), 
c) — > 1 10), and \d) — > |11). Running the Hamiltonian 
H{t) along certain time intervals for appropriate choices 
of T provides a realization of holonomic conditional gates 
of the form 

U(C ,Ci) = |0> (0| ® [7(C ) + |1) (1| ® C/(d). (4) 

Alternatively, one can see that by switching the roles of 
control and target qubits, which can be done by swapping 
the encoding of states \b) and |c), holonomic conditional 
two-qubit gate of the form U(C ) <g> |0) (0| + U(d) <g> |1) (1| 
can be carried out. 

We now verify that C/(Co,Ci) is sufficient for imple- 
menting a universal set of one- and two-qubit gates. 
The key observation is that any triplet Uq,Ui, and D 
can be realized since T is arbitrary up to the condi- 
tion detT ^ 0. First, by choosing T > such that 
cos (a r Z?) = Z, it follows that U\ = Uq and U(Co,Ci) 
reduces to the one-qubit holonomic gate UqZUq span- 
ning the full Lie group of SU(2). Thus, any single- 
qubit transformation is possible utilizing appropriate se- 
quences of T's. Secondly, by choosing T with U\ ^ 
Uq and cos (a T D) = Z implements the holonomic gate 
U(C ,C-l) = |0) (0| ® UqZUI + |1) (1| ® UxZU\ that may 
entangle the qubits. Any entangling two-qubit gate is 
universal when assisted by arbitrary single-qubit trans- 
formations 
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In the remaining part of the paper we discuss three pos- 
sible physical realizations of our scheme in feasible solid- 
state systems. The first example consists of four quantum 
dots arranged in a closed ring geometry. The system that 
we have in mind is the well-known two-dimensional elec- 
tron gas supported, e.g., by a GaAs/AlGaAs heterostruc- 
ture, depleted by a convenient set of top gate electrodes 
to form four quantum dots where electrons are confined. 
Devices of this type can be routinely fabricated nowadays 
with state of the art electron beam lithography [lj|. In 
this system, electrons can tunnel from one dot to an- 
other with a tunnel coupling controlled individually by 
inter-dot gates. We also assume that the local potential, 
and hence the available single-particle energy levels in 
each quantum dot, can be individually manipulated with 
back gates. HQC is realized when this four-quantum-dot 
device contains just a single electron. If only one single- 
particle orbital in each quantum dot is relevant, the sys- 
tem can be described by a four-site tight-binding model, 
forming a closed ring. The second-quantized Hamilto- 
nian reads 



H 



4 

E 

fc=l,<7 



4 

fe=i 



,(5) 



where c kcr (c k(7 ) are electron creation (annihilation) op- 
erators for Wannier orbitals localized at dot-site k, with 
on-site energy e k and spin quantum number a =f, wc 



assume periodic boundary conditions — c\ M . The 
tk.k+i are nearest-neighbor hopping parameters, which 
can be controlled by gate voltages. 

In the presence of an external magnetic field B, 
the Hamiltonian in Eq. ([5]) acquires the Zeeman term 
— gB J2kcr ac la c k(T> where the Bohr magneton has been 
absorbed into g. Furthermore, the magnetic field gener- 
ates a gauge-invariant modification of the tight-binding 
hopping parameters via the Peierls phase factors (Hf 



tk,k+l —> Jk,k+1 — ife,fe+ie 



(6) 



where a k .k+i = fk^k+i A • dl is the integral of the vector 
potential along the hopping path from site k to site k + 1. 

We fix the on-site energies tk = by adjusting the 
voltage on the back gates of the device. We now as- 
sume that hopping parameters can be turned on and 
off with the same time-dependent scaling function f2(f) 
by sweeping the top gates. In the invariant Fock sub- 
space spanned by the four (spin-polarized) one-electron 
states {||000) , 00^0) , |0f00) , |000f)}, the time depen- 
dent Hamiltonian in this ordered basis has the form of 
Eq. ([1]), with the time-independent 2x2 submatrix 



T = 



Jl,2 

J"* 
2,3 



4,1 
J3A 



(7) 



where we have put Jk,k+i — fi£l(t)<Jk,k+i- Note that for a 
given magnetic field we must choose the complex- valued 
hopping parameters so that their phases ak,k+i add up 
to the total magnetic flux $ = § A • dl = AB, where 
A is the area enclosed by the four-dot ring. A change 
of these phases corresponds to a gauge transformation 
A — > A + VA under which the holonomies of the two 
subspaces I = 0, 1 transform into U(Ci) — > TiU(Ci)t1, 
where T ; = diagje"^ A < /+1 ) , e^tMH-s)} with A ( fc ) indi _ 

eating the value of A at site k. This shows the holonomies 
associated with this field-dependent tight-binding Hamil- 
tonian are gauge covariant quantities. 

In order to take advantage of non-adiabaticity to con- 
trol decoherence effects, the typical switching time r 
should be as short as possible. At the same time, the 
action of the Hamiltonian must be " adiabatic" enough to 
prevent transitions to higher orbital levels not included 
in Eq. ([5])- This requirement is expressed by the con- 
dition r s> h/Ae, where Ae is the single-particle level 
spacing. For a quantum dot with a size of 100 nm, typi- 
cally Ae ~ lmeV, which would permit r to be as fast as 
10 ps. This time is considerably shorter than the phase 
coherence time for charge qubits in quantum dots (~ 1 
ns 
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As a second physical example of our scheme, we start 
again with the four-dot system introduced above, but we 
consider the case where four electrons are confined in the 
device. Now, electron-electron interactions must be in- 
cluded and we can do so by adding to the Hamiltonian of 
Eq. ([5]) an onsite Hubbard repulsion term U rik-friki, 



3 



where n ka = c krT c ka-- l n ^ ne limit t k ,k+i/U <C 1 double 
occupancy in each dot is suppressed and the low-energy 
physics of the system can be described by the spin Hamil- 
tonian UM 



H s = ^ Jk,k+lSk ■ Sfc+i + ^ Sfc • Qk,k+lSk+l 
k=l k=l 

4 

+ ^D fe , fe+ i • (s fe X Sfc+i) , (8) 
fc=l 

where s k is a spin— ^ vector operator localized at site 
k. In Eq. @, the first term is an isotropic Heisen- 
berg interaction, with antiferromagnetic exchange cou- 
pling Jk.k+ii the second term is an anisotropic exchange 
described by second rank symmetric tensors &k,k+i-> and 
the last term is an antisymmetric Dzialoshinski-Moriya 
(DM) interaction. The general exchange coupling arises 
from the interplay between the electron-electron repul- 
sion and spin-orbit interaction. For the geometry of four 
coplanar quantum dots in the xy plane only z compo- 
nents of the DM exchange vectors Dfc^+i are relevant. 
In an in-plane electric field, the parameters in Eq. (|SJ) 
can be modified so that the effective Hamiltonian of the 
spin-model is given by 

4 4 

H a = 51 s fc • Jfc ! fc + is /c+1 + ^ D k,k+i ■ ( s fe x s fc+i) z • (9) 



k=l 



fc=l 



Here, we have Jk,k+i = diag{ Jk,k+x> Jk,k+i, 0} with 
J k ,k+i ~ t 2 k : k+i/ u and D Lk+i ~ tk,k+i>^k,k+i/U, where 
Xk,k+i is a spin-dependent hopping term describing the 
spin-orbit interaction. 

Let |t) j. (\i) k ) be the eigenvector of the z compo- 
nent of the spin operator at site k corresponding to 
the eigenvalue | (— 4). The four-dimensional subspace 
with S* otai = 1, spanned by the ordered basis M. = 
lli-ttt) , Ittlt) , Itltt) , |ttt4-)} is an invariant subspace 
[17| . To realize the HQC scheme, the parameters Jk,k+i 
and D k k+1 have to be turned on and off with the same 
time-dependent function tt(t). Under this condition, the 
Hamiltonian in the ordered basis A4 has exactly the form 
of Eq. (TT|) with the time-independent submatrix 



T 



Jl,2 



m,2 

iD. 



2..S 



in 



3,4 



(10) 



where we have put Jk,k+i — hQ{t)Jk.k+i and D 



k,k+l- 



k,k+l ~ 

It is important to stress that since both 
Jk.k+i a n d D k k+1 are proportional to the inter-dot tun- 
neling constants, they can be efficiently manipulated with 
time-dependent gate voltages. 

It is worth mentioning that the spin model described 
by the Hamiltonian in Eq. ([9]) can be realized in 
other settings, such as magnetic impurities on the sur- 
face of topological insulators 18], cold atoms in a 
two-dimensional optical lattice [19j, and nuclear spin 



qubits interacting via a semiconductor-hetero junction 
[20J. These systems may provide alternative platforms 
to implementing our holonomic gates. 

The spin Hamiltonian in Eq. ((8]) also describes the 
magnetic properties of a class of single-molecule magnets 
(SSMs) composed of antiferromagnetic spin rings [Toj . 
Such systems would in principle be more scalable, com- 
pact, and reproducible than the devices based on quan- 
tum dots, since their exchange constants and spin coher- 
ence properties are controlled by the molecular structure 



of the molecule, which can be chemically engineered [21 



The systems that we are interested in are triangular 
SMMs such as Cu3,Vis,Co3 [22]]. Their magnetic core 
consists of three s = 1/2 spins positioned at the vertexes 
of an equilateral triangle and coupled by an antiferromag- 
netic Heisenberg exchange interaction. The ground state 
(GS) manifold of these frustrated SMMs is given by two 
degenerate total spin S = 1/2 doublets (M = ±1/2) of 
opposite spin chirality \ = ±1- The energy gap Aj to the 
first excited S = 3/2 quadruplet is typically in the order 
of 1 meV. A spin-orbit induced DM interaction lifts the 
degeneracy between the two chiral doublets with a split- 
ting Aso — 0.1 A j 23j,|24[. What makes these triangular 
SMMs interesting for quantum manipulation is that, due 
to the lack of inversion symmetry, an electric field in the 
xy-plane of the molecule couples the two GS doublets 
of opposite chirality [3, [2^, [25| ■ In the four-dimensional 
GS manifold spanned by the states \\ = ±1, M = ±1/2), 
the effective low-energy spin Hamiltonian in the presence 
of external electric E and magnetic B fields is [23j 



H, 



off : 



AsoC*S 2 + dE-C|| +B -§S, 



(11) 



where C = ({C,, = C x ,C y },C z ) and S = (S X ,S V ,S Z ) 
are respectively the chirality and spin operators. Due to 
the symmetry of the molecule g = diag{gu,gu, g±}. The 
parameter d is the effective electric dipole coupling, and it 
gives the strength of the coupling between the two states 
with opposite chirality brought about by the electric field. 
In CU3 SMMs d is not small .25J , and for typical electric 
fields generated by scanning tunneling microscope (« 10 2 
kV/cm) the spin-chirality manipulation (Rabi) time is 



10 - 10 d ps [23, 2 



To implement our holonomic scheme in this setting, we 
assume B and E are turned on and off by a common scal- 
ing function Sl(t) and that B z is negligibly small. Under 
these conditions, the Hamiltonian in Eq. (|11|) takes the 
form 



A. 



S( ) 



I 
-I 



(12) 



in the ordered basis {|+1, +1/2), |-1, -1/2), |+1, -1/2), 
— 1, +1/2)} of the ground state subspace. Here, T\x = 
T 2 * 2 = \(B X — iB y ) and T12 = T" 2 *i = ~ ^y)i where we 
have put g\\B k = Ml{t)B k and dE k = h~Q(t)£ kl k = x,y. 
This expression of the Hamiltonian shows that triangu- 
lar SMMs such as CU3 can be used to implement efficient 



4 



and robust gates. Although the restricted form of T asso- 
ciated with this system does not allow for universal com- 
putation, we would still be able to obtain some of the 
most important gates such as all the Pauli-matrices and 
the controlled-Z gate in a fashion that is highly resistant 
to errors. The crucial challenge relies on the ability of 
generating electric and magnetic fields characterized by 
time-dependence function Q(t). Fig. Q] shows the fidelity 
F of the Pauli-Y gate performed with a triangular SMM, 
as a function of the spin-orbit induced splitting Aso and 
operation time r. Here, we have adopted the definition 



F = F{W, V) 



|ReTr (WW) , where W is the ideal 



gate and V is the gate carried out by the perturbativc 
Hamiltonian in Eq. §T2$. The fidelity F(W,V) deter- 
mines the Hilbert-Schmidt distance between W and V. 
For typical values of Aso ~ 0.02 meV Q the fidelity is 
over 98% for r = 10 - 100 ps. 



using non-adiabatic quantum holonomies [5| acting on 
qubit pairs encoded in a four level structure. We make 
use of the two-qubit system as a unit cell for scalable 
all-geometric quantum computation. We propose three 
different realizations of our scheme based on an uncon- 
ventional use of quantum dot and single-molecule magnet 
devices. We emphasize that the physical realizations of 
HQC presented here differ substantially from previous 
proposals of implementation of charge and spin qubits in 
quantum dots and single-molecule magnets. Finally, we 
would like to point out that the four level configuration 
proposed here is quite a general scheme that can be real- 
ized with other physical systems actively considered for 
quantum computer implementation T3- 20|. 

This work was supported by School of Computer Sci- 
ence, Physics and Mathematics at Linnaeus University 
(Sweden). E.S. thanks National Research Foundation 
and the Ministry of Education (Singapore) for support. 
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FIG. 1. Fidelity F of the holonomic Y-gate for a triangular 
SMM as a function of the spin-orbit induced splitting Aso and 
operation time r. Here we have chosen g\\ (B x —iB y ) = 
d(E x — iE y ) — —~Kl(t), where Q(t) is a square pulse such 
that f Q T 9,(t)dt = 7T. 

The four-level configuration presented above can be 
used in a scalable pattern necessary for universal quan- 
tum computation on any number of qubits. This can be 
achieved by implementing the Hamiltonian 



7 1 (p</) 

J"(P9)t 



(13) 



where p and q label the qubits. Any gate acting on qubits 
p and q can be implemented by turning on and off tt pq (t) 
while the other terms are kept off. H(t) can be realized 
with the first two examples using the techniques intro- 
duced in Ref. [2(| to design scalable architectures. This 
Hamiltonian would also be realized by arranging trian- 
gular SMMs on a surface in a way that each of them can 
be manipulated individually. 

In conclusion, we have introduced a scheme to imple- 
ment high-speed universal holonomic quantum gates that 
considerably differs from the other setups proposed for 
HQC. It allows for implementation of conditional gates 
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